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^ . Abstract 

D ■ We calculate the two-pion exchange potential of DD and DD by using Gaussian type form factor to regularize 

' loop integrals. Then we discuss the possibility of the molecule states composed of D mesons. Our results show 

it is very difficult to bind DD by pure pions exchange. On the contrary, the possibility of the DD molecule 
state existence cannot be completely ruled out. We also extend our discussion to BB and BB . 
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^ ■ 1 Introduction 

^ . The idea of molecular state of mesons was first raised in [Tl[2], and it was systematicly studied by N.A.Tornqvist 
QJ ' in [3]. This idea is often used to study the resonances with "strange " decay property as well as a mass closed 
,. to a respective two-meson threshold. One of the newest example is the hadron state of X(3872). A lot of papers 

suggested it to be a molecular state of D'^D* . 

The X(3872) was first discovered by Belle [4j in J/i/jTr+Tr" channel, and soon confirmed by CDF [5 , DO [6 ,and 
^ ' BaBar f^. Its mass is Mx = 3871.2 ± 0.5MeV and width is F < 2.3MeV. The study of the angular distributions 
' and correlations of the J/-07r+7r~ final state shows that J^'~^ — 1++ and J^'-' — 2 are two possible quantum 
' number assignments. Since the 2"^ Pi cc state x'ci in quark model is at least 50MeV heavier than X(3872) and the 
■ decay of a pure charmonium state into J/ipiriT violates isospin symmetry, it is very difficult to assign X(3872) into 
[ charmonium. Many new reasonable explanations were raised, such as explaining X(3872) as a molecule state , a 
hybrid charmonium , a diquark anti-diquark bound state, a tetraquark state or so. Since X(3872)'s mass is just a 
' half MeV below DD* mass threshold and its branch width in the DD* channel is ten times larger than that in the 
J/^Tr~^Tr~ channel, the molecule picture becomes one of most favorable model though this picture meets a problem 
in explanation of the decay X{3872) — > J/V'tt+tt^. 
, . The molecule picture for X(3872) have been discussed in many papers [51 [HI UHl [HI [ID- As the proximity of 

' bound state is always built in a potential model. In [10] Swanson suggested the potential is mainly supplied by 
exchanging pion mesons and quarks since the one pion exchange is not good enough to bind D'^D together. In [13j 
the authors traced back to pion exchange potential of deuteron in nuclear physics, where the single pion exchange 
cannot bind proton and neutron either. Only by introducing the sigma meson exchange would the deuteron to 
be stable. So they had introduced the sigma meson exchange in the D^D . However their result showed the 
situation became worse. A bound state of D^D* could appear only when the parameter of energy cutoff A become 
unphysically large. Their results disfavors the X(3872) being a bound state of DD*. In [14] C.E.Thomas and 
F.E. Close calculated the pion exchange of D and D*m the quark level. And they did find a bound state in their 
potential. And the bind energy is sensitive to the energy cutoff parameter A and the coupling constant g. Beside 
the pion exchange, other mesons exchange, such as cr, p,uj and 77 exchange between D and D* , has been discussed in 
[T5] . The authors found a loosely bound state of DD* when the energy cutoff A = llbOMeV. Their results favors 
the X(3872) being a loosely bound state of D and D* mesons. 

Until now, from both of the theory and the experiment, there are two opposite opinions about the molecule 
picture of the resonance X(3872). However, more new resonances, such as X,Y, Z series of the charmonium, are 
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found continuously in experiments. Some of them are also difficult to be explained as the pure cc states. Therefore, 
regardless of whether the X(3872) is a molecule state of DD*, the interaction between D mesons in hadron level 
will be still interesting. 

Our work starts from the thinking of the two-pion exchange potential between D and(or) D* mesons. In nuclear 
physics, the two-pion exchange potential plays an important role in explaining the inside structure of nucleons. 
We guess such effects are also important between the D and(or) D* mesons though it was argued that the chiral 
perturbative expansion converges more quickly in heavy mesons. At first, we calculate the two-pion exchange 
potential between D and D by using Gaussian type form factor to regularize loop integral and study the bound 
state of DD. We find, though there is not the one-pion exchange potential, a stable bound state of DD can exist 
at A = 1750MeV. Then we perform the similar calculation for the two-pion exchange potential between D and 
D* . Meanwhile, we also use the same form factor to regularize the one-pion exchange potential between D and D* 
for comparison. We find the chiral perturbative expansion is quite convergent even when the cutoff A ~ 5GeV. 
This result shows the regularization scheme works quite well. Similarly, we use the obtained potential to study the 
bound state of DD* and find a stable bound state exists only when A > A.JGeV. This result is better than [T5] 
though it is still negative for the existence of the molecule X(3872). Our results show, even in heavy mesons, the 
multi-pion exchange is still important. 

We organize our paper as the following. In section 2, we give the effective Lagrangian that we use in the 
calculation. This Lagrangian is derived from the most general one given by Mark.B.Wise in [TB]. In section 3, we 
first calculate the two-pion exchange potential between DD, then we numerically solve the schrodinger equation to 
find whether this potential can form a bound state. In section 4, we calculate the potential of DD and discuss the 
possibility of the existence of the molecule X(3872). In section 5, we extend our discussion to the case of the B and 
B mesons . Finally we give a brief summary and conclusion. 

2 Effective Chiral Lagrangian 

The effective chiral lagrangian for heavy mesons was first given by Mark . B . Wise 1 1 6j : 

£ - -iTrHai^^d'^Ha + h.i:THaHbV^^{£}d^i + m'^)ba + UgTvHaHbiul^{£}d^^ - id£})ba + • • • , (1) 
where 

Ha=^-^{P:^r -Pal^) (2) 

and the P and P* are meson fields and ^ is the nonlinear realization of the SU(3) goldstone fields. 
For our purpose, we can expand the lagrangian JT]) as 

C = -^Pt-{TTXTr)P^j-{P*^V2t-^TTP + P^P*V2t-^n)-^P*t-{Trxn)P*, (4) 

where t — ^a, f = ISOMeV is the tt constant, g is the coupling constant(we use g — 0.59 in the calculation 
below). Same as in [TCj, the factors of V^Moand \/Md-' have been absorbed into the P and P* fields, and the 
Lagrangian is independent of the heavy-quark mass. 

3 DD Potential and Bound State 

Now we start to calculate the potential of D and D. Since there is no one-pion exchange between D and D mesons, 
we are directly going to calculate the two-pion exchange potential. There are four two pion-exchange Feynman 
Diagrams as showed in Fig.l. 
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Fig.l The parallel box(a), crossing box(b), triangle(c) and two-pion loop(d) digram. 

Our calculation is very similar to that for the nucleon-nucleon potential [T7] [3T], except that the parallel box 
diagram in our case is a little bit different. In the nucleon case, it should be careful to extract " iterated one- 
pion exchange" contribution 18j. Therefore, the "old-fashioned" time-ordered perturbation theory is widely used 
[T7] . However, there is not one-pion exchange between D and D mesons, we can calculate directly using covariant 
perturbation theory. 

The mass difference A — Md* — AId plays an important role in our calculation, so we write the propagator of 
the D* meson like ( for example of Fig. la) 



-k - 2A + ie' 



(5) 



where we use the conditions pf ~ and qo = 0. Then the rest calculation can be carried out directly like in [T7] . 
3.1 Two-pion Exchange Potential 

As usual, we work in center-of-mass(c.m.) system. And denote p as the intial momentum, p' as the final momentum, 
q = p — p' as the transferred momentum. In the following, we will use the symbols: 
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a = 2A = 2{Md- - Md) > 

After straightforward calculating each figure in Fig.l, we can get the potentials in momentum space: 
for the diagram a in Fig.l, parallel box digram: 
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1 ^ff^4, 
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0^1(^1 + 0) a;2('^2 + a) 
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the diagram b in Fig.l, crossing box digram: 
the diagram c in Fig.l: 
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1 



the diagram d in Fig.l: 



1 25 2 2 f dH 

16 P^''^^J (2^ 



1 / 1 



1 



2 \a;i -|- a W2 + a/ — 
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(8) 



(9) 
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where ti and t2 are the isospin matrices in the line 1 and 2 respectively. 

To transform the potential from momentum space into coordinate space, the momentum space cutoff is used. 
This is a conceptually simple way and is mathematically easy too, we just need to bring a cutoff function and 
choose a suitable cutoff energy. Although it is not theoretically exactly true in physics, but this has been used 
widely in phenomenology physics, and is believed to be consistent with how the real world is going on. Like in 
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the reference [TOl [20] HH, we add a Gaussian form factor exp{—-^) in each vertex (in momentum space ) in the 
diagrams of Fig.l, that can regulate the loop integrals in the potential. Thus, the integrals can easily be carried 
out numerically. When we do the transformation, we use the formulas and techniques in Ref [21]. Here we only 
give the simple potential formulas. 

After the Fourier transformation, the potential between DD can be obtained as following( where the function 
F(X,r) can be found in Appendix): 

V{r) = Va{r) + Vb{r) + K(r) + Ki(r), (13) 

where 



Va{r) 
Vb{r) 
Vc{r) 



dX 



~i:-il)H3-Sh-h) [ dX- 

iDTT / 

l.l.t,.t,j dXX^FWr) 



a2 + A2 



A2 



(a2 + A2)2 



\f' (A, r) f' (A, r) + f" (A, r) f" (A, r) 



\f' (A, r) F (A, r) + f" (A, r) f" (A, r) 



A2 



F {X,r)F (A,r) 



(14) 
(15) 
(16) 
(17) 



The potential is a function of energy cutoff parameter A. In Fig. 2 we show the potential with different A. We 
can see the potential is sensitive to the cutoff A. The value of A is not known right now, however, it is generally 
believed that A is around IGeV. At most, it should not be larger than the mass of the D meson. We will discuss 
its value in the following. 

It should be also mentioned here that, what we deal with the DD potential is also available for the DD potential, 
i.e., there is no difference between DD and DD up to the order of the two-pion exchange . Therefore, in the following 
what we talk about DD is also valid for DD. 
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Fig. 2 The potentials of DD with various A. The left is in the channel 1=0 and the right is in the channel 1=1. 
3.2 Numerical Results 

We numerically solve the Schrodinger equation to find the bound state of DD. Our result is certainly dependent 
on A. For the channel 1=0, we find that a stable bound state exists only when A > ll^QMeV . For 1=1, a stable 
bound appear after A increases to 2750MeV. We summarize our results in Table. 1. Although the potential in both 
of the channel 1=0,1 is attractive, it seems that it is not deep enough to bound a molecule state if we believe A 
should not be larger than IGeV. However, the experimental data about nucleus support a larger A = llSOAfeF 
[20] • Considering the mass of the D meson is heavier than nucleon's, it could not completely exclude the possibility 
of the existence of the DD molecule state. 



A[MeV] 


Eund [Mev] 


E = M - Eund[MeY] 


A[McV] 


Ebind [Mev] 


E ^ M - Eb.,nd[M^eY)] 


1700 


Not found 




2700 


Not found 




1750 


-4.54 


3725.12 


2750 


-1.53 


3728.13 


1800 


-21.1 


3708.56 


2800 


-14.71 


3714.95 


1850 


-48.57 


3681.09 


2850 


-36.19 


3693.47 


1900 


-90.68 


3638.98 


2900 


-68.68 


3660.98 


1950 


-151.04 


3578.62 


2950 


-114.53 


3615.13 


2000 


-233.54 


3496.12 


3000 


-176.06 


3553.60 



Table. 1 The binding energy and the corresponding A for the DD bound state. The left one is in the channel 1=0, 
the right one is 1=1. 



4 DD* Potential and Bound State 

The calculation of the potential between D and D up to the order of the two-pion exchange can be carried out 
in the same way in the framework of the lagrangian The only exception is that the one-pion exchange plays 
an important role in the case of DD* . The existence of the one-pion exchange also makes the treatment of the 
parallel box diagram different. After giving the potential, we will discuss the possibility of the existence of the D D 
molecule state. 

4.1 Two-pion Exchange Potential 

Since there is the one-pion exchange between D D , it needs to extract the "iterated one-pion exchange" part from 
the parallel box diagram. Normally, such extraction can be done by using the time-ordered perturbation theory. A 



more convenient way for us is done in the covariant perturbation theory. As pointed out in [18j . in the covariant 
formuhsm , the "iterated one-pion exchange" part in the paraUel box diagram can be extracted by simply omitting 
two poles, which is moving towards each other in heavy mass limit, in the /q integration. Then we obtain the 
two-pion exchange potential of DD in the momentum space: 



VaD'D 
VbD'D 

VclD'D 
Vc2D'D 
VdD'D 



167' 



(3-8ti -ta) 



J37 1 



1 



1 



1 



j3; 1 



1 



K-a2)(^2_„2) 



1 



a;i(wi+a) W2(w2 + a) 



1 



= 7(7) 72*1 • t2laP 



(2^) 



2a 



47' P 



1 1 

4? 



dH 
(2^ 



a;| — oj^ \^ cji + a a;2 + a 



W2 



w| — wi + a a;2 + a ^ 



(18) 



The transformation to coordinate space can be done in the same way as we did to DD. Then we can get the 
two-pion exchange potential in coordinate space: 



VDD'{r) = VaD'D{r) + VbD'oir) + VaclD'oir) + Vc2D'D{r) + VdD'D{r), 



(19) 



where 



VaD'oir) 
VbD'D{r) 



ID'D 



Vc2D'D{r) 

VdD'oir) 



(a2 + A2)2 



;F iX,r)F (A,r)+F {X,r)F (A,r) 



1 1(^)4(3 _8ii.t2)(/- 
16 TT / 

-Y^(fr(3 + 8i,.t2)(/-^:^)„,- 



2a 

TT 



dX 



a2 + A2 



;F (A,r)F (A,r)+F (A, r) F (A,r) 



} 



i^rlh-t2{i-'^^ 



i(ff)2. 
47^ /2 



)c/3-[2ae^/2(V4TO2;^7i2, r)^ ( Virnf^T?, r) 



/ dX- 



A2 



F (A,r)F (A,r) 
A2 



a2 + A2 



\{^j?j,h-t2l^plj 

~ti-t2lap- f dXX^F\X,r), 

^ I J 



^' (A,r) f' (A,r)) 



(20) 



where s is the spin matrix of the D* and a/3 is the indices of s. 

Since a — 2 A > 2m7r, the term a/ 4m2 — a2 in above formulism is a imaginary number. This makes the potential 
()19p to be complex . In fact , the one-pion exchange potential of DD* is also complex. In the following we only 
choose the real part of Vdd' (r) in solving the Schrodinger equation as many authors did in the studies of the DD* 



molecule state. Besides, we will average the matrix (/ - 



'-)ai3 in the numerically calculation. 
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4.2 One-pion Exchange Potential 
The one-pion exchange potential is: 

ViAq) = - J(7)'2ti • ij-— (21) 

3 / + — A'' 

In the transformation to coordinate space, we also use the Gaussian type form factor for consistency . The 
corresponding potential is 

V^Ar) = -lijf^t^ ■ hA'^Wml - A2,r) (22) 

Since the monopole type form factor F{q) = \f^i is used widely to regularize the one-pion exchange potential 

of DD , it is deserved to compare these two different type form factors by varying A. Wc find A = lOOOMey in the 
Gaussian type is about A = 900MeV in the monopole type, therefore there is almost no difference between these 
two types form factors in the region of A where we are interested in. Before proceeding the discussion of the bound 
state, it is also useful to consider the convergence of the chiral expansion. In Fig. 5, we show the dependence of the 
one-pion and the two-pion exchange potentials on the cutoff A. In fact, even when A = 5GeV, the chiral expansion 
is convergent quite well. So the regularization scheme works well. 

4.3 Bound State 



A[McV] 4650 


4700 


4750 


4800 


4850 


4900 


4950 


5000 


Ebi„d[MeV] Not Found 


-0.99 


-2.58 


-4.61 


-7.14 


-10.23 


-13.91 


-17.90 



Table. 2 The binding energy of the DD* bound state in the channel 1=1 and the corresponding A. 

In Table. 2 we summarize the bind energy of the DD* bound state and the corresponding value of A in 1=1 
state. One can find a stable bound state appears only when A > 4700MeV^. So the pure pion-exchange cannot 
bind D and D* as the X(3872). 

5 B and B* Bound State 

It is easy to extend above results to the case of B and B* . The only difference is that the mass difference between 
B and B* is smaller than the mass of tt, so all potentials are real. The values of the B and B* mesons in our 

calculation are Mb = 5279. 5Mey and Mb' = 5325.17^61^. These masses are much bigger than the masses of D 
and D*, the bound states of BB and BB are more possible to exist. We list our results in Table. 3 and Table. 4 



A[McV] 


Ebind [Mev] 


E = M - EbindlMeV)] 


A[MeV] 


Ebind [Mev] 


E = M - EbindlMeV)] 


1450 


Not Found 




1600 


Not Found 




1500 


-1.59 


10557.41 


1650 


-1.71 


10557.29 


1550 


-16.71 


10542.29 


1700 


-6.72 


10552.28 


1600 


-50.42 


10508.58 


1750 


-16.06 


10542.94 


1650 


-107.32 


10451.68 


1800 


-31.03 


10527.97 


1700 


-193.49 


10365.51 


1850 


-53.09 


10505.91 



Table. 3 The bind energy of the BB bound state and the corresponding A. The left one is the 1=0 state and the 
right one is the 1=1 state. 
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A[MeVJ 


Ebind [Mev] 


E = M - Ei,ind[MeY)\ 


1700 


Not Found 




1800 


-0.12 


10604.58 


1900 


-1.47 


10603.13 


2000 


-3.98 


10600.62 


2100 


-7.88 


10596.72 


2200 


-13.36 


10591.24 



Table. 4 The bind energy of the BB* bound state in the channel 1=1 and the corresponding A. 

6 Summary and Conclusion 

In this work, we calculate the two-pion exchange potential between D mesons and(or) D* mesons by using Gaussian 
type form factor to regularize loop integrals. We check the convergence of the chiral perturbative expansion by 
varying the cutoff A and make sure the validness of the regularization scheme. By using the obtained potentials we 
discuss the possibility of the existence of the heavy-meson molecule states. Opposite to the naive expectation, we 
find the existence of the DD molecule state is more possible than that of the DD* molecule state. The existence 
of a DD molecule state requires a large cutoff A = 1750MeT^. Comparing with A = llSOMeT^ in the deuteron 
case[20 , a DD molecule state indeed is not favorable. However, if we consider the fact that the D meson is heavier 
than nucleons, we cannot rule out completely the possibility of the existence of the DD molecule states. Similar 
to [13], our result excludes that the pure-pion exchange can bind D and D* mesons as a molecule state. We also 
extend our study to the case of B and B* mesons. Since the mass of B meson is much heavier than the D meson's, 
the existence of the B and(or) B* molecule state is more possible than that in the D meson case. However, the 
situation does not become much better than that of the DD case. 
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Appendix A 

Here we list the functions used in the potentials: 



d_ 

dr 

dr 



F{X,r 
him, A 



e^/2(V(2nvF+A2, r) 
m 
47r 



(t)°{m,r) 



e 



Ar m 



. , e"'''Erfc( 
A' ■' ^ 2 



Ar m . 



A' 



2mr 



A r ni A V 777 

et^ <; (1 + mr)e-'"''Erfc{ — - + ^) " (1 " mr)e"''' Er f c{— + —] 



2 



2(mr) 



= 0°(?77,,r) 



^/A 



—m r(f)!.g{m, r) 

m^4>l(m,r) - --^0°(m,r), 
r dr 



where the function Erfc{x) is the complementary error function. More details can be found in the Appendix 
of Ref 211. 



Appendix B 



Here we show the figures of two-pion and single-pion exchange potential of DD 
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Fig. 3 In this figure, we show the potential of DD in the channel 1=1 when A = lOOOMey. The real line is the 
two-pion exchange potential, the dashed line is the one-pion exchange potential, and the dotted line is the total 
potential. 
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Fig. 4 In this figure, we show the total potential of DD in the channel 1=1 with different A. 
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Single-pion exchange potential 
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Fig. 5 Here we show the single-pion exchange potential and two-pion exchange potential of DD in the channel 1=1 
with different A. 
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